Comment on "Exact results for survival probability in the 
multistate Landau-Zener model" 

B. E. Dobrescu 
Department of Physics, Texas A&M University, 
College Station, TX 77843-4242, USA 

N. A. Sinitsyn 

Department of Physics, University of Texas at Austin, Austin TX 78712-1081, USA 

Abstract 

We correct the proof of Brundobler-Elser formula (BEF) provided in [2004 J. Phys. B: At. Mol. 
Opt. Phys. 37 4069] and continued in Appendix of [2005 J. Phys. B: At. Mol. Opt. Phys. 38 
907]. After showing that some changes of variables employed in these articles are used erroneously, 
we propose an alternative change of variables which solves the problem. In our proof, we reveal the 
connection between the BEF for a general iV-level Landau-Zener system and the exactly solvable 
bow-tie model. The special importance of the diabatic levels with maximum/minimum slope is 
emphasized throughout. 
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I. INTRODUCTION 

In Ref. jjj Volkov and Ostrovsky (referred to below as V-O) proposed a rigorous mathe- 
matical proof of the Brundobler-Elser formula (BEF) [3( for a general iV-level Landau-Zener 
(LZ) system. The approach of |l| was to consider a formal solution in terms of an infi- 
nite series in powers of the coupling constants, with further calculation of every term and 
summation of the resulting algebraic series. Following this strategy, these authors have 
demonstrated the exact cancellations of some quantum paths and have reduced the problem 
to the calculation of a sub-series (5.1) (the numbering of formulas is as in Q]). 

In the calculation of this remaining sub-series (5.1) in |l| V-0 made the error of assuming 
a wrong domain of integration for their variables Xj (see Eq. (5.2) of |l|). 



In Ref. |2| , whose main text is dedicated to a different subject, V-0 included an Appendix 
in which a new attempt is made at calculating the sub-series (5.1). After acknowledging 
the error made in l|, V-0 proceeded to employ a different change of variables, previously 
proposed by Kayanuma and Fukuchi in [5j. 

We would like to point out that this derivation contained in the Appendix of Q] is 
still erroneous. After Eq. (A. 4) (in the notation of Ref. Q|) the authors claim that the 
integrand is symmetrical under the permutation of variables xo, ■ ■ ■ , It is, however, 

straightforward to check, by direct examination of the terms containing the 4-th power in 
the coupling constants Vi, that this is not true; nor is it symmetrical under exchange of pairs 
of variables (x , yo),---, (x m -i, Um-i)- 

The following example will support these considerations: the term corresponding to I = 



m = 2 in Eq. (5.1) of 



reads 



e i^ii 2 e i^ 2 i 2 r dt ± r r ^ r ^ ^ 
fe^i k 2J ti 

x exp [i (si - e kl ) (U -t 3 )+i (ei - e k2 ) (t 2 - h)] 
x exp [| (A - kl ) (t\ -tj) + % - (A - (3 k2 ) (t\ - tfj 

Let us consider the simplest model beyond the two-level LZ system, namely the one with 
N = 3. Then, the sum above has 4 terms 
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dt 2 / d£i 

-oo J — oo 



x exp [i (ex - e 3 ) (t 4 - t 3 ) + i (s, - e 2 ) (h - h) 
~\ (A - A) (tf - tf) + \ (A - A) (3 - *?) 



x exp 



The change of variables p used by V-0 in Appendix of |2j amounts to 



(5) 



ti = Xi, Xi G (— oo, oo) 

t 2 = X! + yi G [0, oo) 

h = x 2 + yi, x 2 G [xi, oo) 

U = x 2 + yi + y 2 , y 2 G [0, oo) 

The integrals I 22 and J 33 can be easily evaluated (see p 1 ), e.g. 



12 / «^2 / dxi / dy 2 / dy\ 

oo J —oo JO JO 



I 22 = \Vi 2 \ I d.r 
x exp 
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and, similarly, J 33 



I Vis I' 
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The sum of the remaining two integrals reads 

/oo roc roc roc 

dxi / dx 2 / dy 2 / dy x (7) 
-oo Jxi JO JO 

x {exp [i (e 1 - e 2 ) y 2 + i (ei - £3) Vi] 

x exp I (J3i - p 2 ) (2x 2 y 2 + y 2 2 ) + % - (pi - fo) (2x lVl + yf) +i((3 l - p 2 ) Vl y 2 
+ exp [i (ex - e 3 ) y 2 + i (si - e 2 ) 

x exp i (ft - /3 3 ) (2x 2?/2 + y 2 2 ) + I (A - /3 2 ) (2x iyi + yfj +i(J3i- /3 3 ) yiya 

Contrary to the claims made by V-0 after Eq. (A. 4) of 0, the integrand of (0) is clearly 
not symmetric with respect to interchanging the pairs (xi,yi) and (x 2 ,y 2 ). Indeed, the 
terms {(3\ — (3 2 ) y\y 2 7^ — /3 3 ) y\y 2 belong to different terms of the sum in the integrand 
and are left unchanged by the permutation. Therefore, one cannot extend the domain of 



integration of variable x 2 to (—00, +00) as proposed by V-0 in Eq. (A. 5) of 



While the mathematical errors made by V-0 are presented above in detail only for com- 
pleteness, there is also a more general argument as to why the calculation of the series (5.1) 

n n 

in both Ref. m and Ref. U is erroneous. One can notice that this series represents a formal 
solution for a special case of the multi-state LZ-model, in which all states interact with only 
one special level with slope Pi. Series (5.1) is the probability amplitude for the system to 
remain on this special level after all crossings if only this level has been initially populated. 
This observation does not assume at all that the slope Pi is an extremum. Therefore, if 
the proof provided by V-0 in [1] and |2j were correct, then their summation of the series 
(5.1) that does not use the assumption of a maximal/minimal slope for the level 1 could be 
applied to such a LZ-model with an arbitrary slope Pi. 

However, it is well known from numerical simulations [3] of this model that if pi is not an 
extremum, then the BEF does not hold, and the transition probability for 1 — > 1 depends 
essentially on the parameters e^. Also, the exactly solvable case of the above example in 
which all £j = 0, the bow-tie model proves that the diagonal elements of the S matrix 
are different from BEF except for the levels with maximal /minimal slope. 

The total ignorance of this very important property of the slope Pi being an extremum 
in the calculation of series (5.1) by V-0 has actually been the trigger that prompted our 
search for their mathematical mistakes. 
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II. PROOF OF BEF 



We will split the proof of BEF in three steps. The first step has been correctly performed 
in and we only briefly mention here the main results in order to introduce the notation 
and starting point for the next steps. 

It has been shown in lj that if the slope (3\ is maximal/minimal then the transition 
amplitude for the transition 1 — > 1 has the form: 

oo N N N 

S u = <l|[/(oc,-oc)|l> = l+£(-ir£ £ \Vi kl \ 2 \V lk2 \ 2 ---\V lk J 2 x 

m=l fci=2fc2=2 k m =2 

/OO l-Tl rT2m-l 

dn / dr 2 ■ ■ ■ / dr 2m 
-oo J —oo J —oo 

x exp [i^M ( rf _ T =) + ( T | _ ^ + . . . + ,*M ( tL _, _ tL ; 

x exp [i£ (fci) (ti - r 2 ) + i£ (fc 2 ) (r 3 - r 4 ) H h (fc m ) (r 2m -i - r 2m )] , 



where 



5 (kj) =fa- (3 kj and J5 (%) = £i - j = 1, 2, . . . , m. (9) 



U is the time-evolution operator and N is the number of levels in the system. 

Series (JEJ) is equivalent to the series (5.1) in [1]. From this point on, Volkov and Ostrovsky 
do not use the property of /3\ being maximal/minimal anymore in |l| and in the Appendix 
of P|, and the changes of variables they employ are used erroneously, as shown in the 
Introduction. 

In the second step we prove that if (3% is maximal/minimal, then Su = (1\U (oo, — oo) |1) 
does not depend on E for any lj = 2,3, N and any j = 1, 2, . . . , m. 
From Eq.fjHJ) it follows that: 

oo N N N 

^ = E(-i) m E E ••• E \Vi kl \ 2 \v lk2 \ 2 ---\v lk jx 

OCj \ L j) m=l k 1= 2k 2 =2 k m =2 

/oo rri rTim-i 

dn / dr 2 - ■ ■ dT 2m 
-oo J — oo J — oo 

[5 (Z j; fci) (n - t 2 ) + 5 ((,-, fc 2 ) (r 3 - r 4 ) H h <5 (/ i; fc m ) (r 2m _i - r 2m )] 

x exp (xf - r|) + ^ (r| - r|) + • • • + (r| m _, - 4,,) 

x exp [zE (fci) (n - r 2 ) + iE (k 2 ) (r 3 - r 4 ) H h z£ (fc TO ) (r 2m _i - r 2m )] , (10) 



where 5 k p ) is the Kronecker delta. 
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Next, we introduce the well-known change of variables (see, e.g., textbook |7[, page 114): 



T\ = X\ G ( — oo, oo) 
Tj + i = Tj — Xj+i, with Xj + \ G [0, oo) , j = 1,2,..., 2m — 1. 

From Eqs. (HUJ) and (HU) it follows that 



N N 



N 



= £(-i) m £ £••• £ I^J 2 l^ 2 l 2 ---I^J 2 x 

<7 - C/ I'jV m=l 



fe 1= 2 fc 2 =2 fc m =2 
/■oo /*oo 

x / dxi / dx 2 \ dx 3 ■ 

/-oo JO JO 



dx 



2)ii 



xi [5 h) x 2 + 5 k 2 ) x 4 H h 5 (lj, k m ) x 2m ] 

xF(x 1 ,x 2 ,x 3 , ...,x 2m ) , 



(11) 



(12) 



where 



F(xi,x 2 ,x 3 ,...,X2m) = exp 



(^2, X4, . . . , X2m-2, ^2r 



x exp [i {B (ki) x 2 + B (k 2 ) x 4 H \- B (k m ) x 2m } x x ] 

x exp [-i {B (k 2 ) x 4 + B (k 3 ) x 6 H h -B (A; m ) x 2m } x 3 ] 

x exp [-z {5 (k 3 ) x 6 + B (k A ) x$ H h -B (& m ) x 2m } x 5 ] 



x exp [-i {B (k m ) x 2m } x 2 



m— lj > 



(13) 



and 



G (x 2 , x±, . . . , x 2m - 2 , x 2r 



B (h) (x 2 2 ) + B (k 2 ) (2x4X2 + 
+B (k 3 ) (2x 6 x 2 + 2x 6 x 4 + 5 



+73 (& m ) (2x 2m £2 + 2x 2m x 4 H h 2x 2m x 2m - 2 + x\ m/ 

-2 [E (h) x 2 + E (k 2 ) x 4 + --- + E (k m ) x 2m ] . (14) 



Upon integrating over x\ in Eq. (jl2j) one obtains: 



OS 



11 



N N 



N 



dE (I,) 



E(-i) m E £••• E l^ri^ 



ifc 2 1 



l^i fem | 2 x 



m=l fcl=2 &2=2 fcm=2 

/•OO /"OO /'OO 

x / dx 2 / (fX3 • • • / (ix2 m 

JO JO JO 



xi [5 ki) x 2 + 5 k 2 )x 4 -\ h 6 (lj, k m ) x 2m ] 



x exp 



— — G (x 2 , X4, . . . , £ 2m _ 2 , X 2r 



x 2%5 (B (kt) x 2 + B (k 2 ) x 4 + --- + B (k m ) x 2m ) 

x exp [-i {B (k 2 ) x 4 + B (k 3 ) x 6 H h B (k m ) x 2m } x 3 ) 

x exp [-i {B (k 3 ) x 6 + B (fc 4 ) x 8 H h -B (fc m ) x 2m } x 5 ] 



x exp [-z (fc m ) x 2m } x 2 



m— 1 



(15) 



From Eq.® it follows that, for any kj, B (kj) > if (3\ is maximal, and B (kj) < if 
/3i is minimal. Therefore, 5 (B (ki) x 2 + B (k 2 ) x± + ■ ■ ■ + B (k m ) x 2m ) = unless x 2 — X4. — 
■ • ' — x 2m — 0. The presence of the term 



[5 (lj, hi) x 2 + 5 (lj, k 2 )x A -\ h 5 (lj, km) x 2r< 



(16) 



in the integrands of Eq. (|15|) ensures that each of the integrals is zero. 

This argument can be made rigorous by regularizing the behavior of integrals at ±00, as 
follows 



/oo /*oo /*oo roo 

dx x / dx 2 / dx 3 --- dx 2m F (x u . . . , x 2m ) x 2q 
-00 Jo Jo Jo 

f roo roo roo roo roo roo 

= lim <^ / dx 2 / dx 4 ■ ■ ■ dx 2m / dx\ \ dx 3 ■ ■ ■ dx 

»?->o+ t Jo Jo Jo Jo Jo Jo 

xx 2g F (x Xl . . . , x 2m ) e -v(*i+x 3 +-+*2 m -i) 

roo roo roo rO roo roo 

+ / dx 2 I dx4 ■ ■ ■ dx 2m I dx\ \ dx 3 ■ ■ ■ dx 

Jo Jo Jo J-oo Jo Jo 

x x 2q F(x 1 ,...,x 2m )e-^ +x:i+ - +X2 -^ }, 



2m- 1 x 



'2m-l x 



(17) 



for any q — 1, 2, . . . ,m. Upon integrating over the Xj with odd j in Eq. (|T7jl . one obtains 

roc roo roo 

I q = lim / dx 2 / dx A ---l dx 2m W (x 2 , x 4 , . . . , x 2m ; rf) , (18) 
v^o+ Jo Jo Jo 

where 



W (x 2 ,x 4 , . . . ,x 2m ;r]) = exp 



-G(x 2 , x 4 , . . . , x 2r 



2r]x 2q 



[B (kt) x 2 + B (k 2 ) x 4 



B(k m )x 2m ] +r] 2 



x 



i [B (k 2 ) x A + B (k 3 ) x s -\ V B (k m ) x 2m ] + r) 
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X i[B (k m ) x 2m ] +rj' ^ 
Since either B (kj) > or B (kj) < for any kj, from Eqs. (fTH|) and (jP3|) it follows that 

roc roc />00 

\L\ < lim / dx 2 / dx4--- dx 2m \W (x 2 , x±, . . . , x 2m ; t])\ 
■n->o+ Jo Jo Jo 

rv/\B{k x )\ p/\B(k 2 )\ p/\B(k m )\ 

= lim / dx 2 / dx^ ■ ■ ■ I dx 2m x 

»7->-0+ Jo Jo Jo 



2r]x 2q 



[B (hi) x 2 + B (k 2 ) x 4 + --- + B {k m ) x 2m ] z + v 2 

1 1 



,(20) 



[B (k 2 ) x A + B (k 3 ) xq H h -B (fc m ) x 2m f + ?? 2 (fc m ) £2m] 2 + 7/ 2 

For any continuous function / defined on a compact interval [a, 6] the relation 
f a f (x) dx = f (£) (b — a) holds, for some £ £ [a, b\. Hence, Eq.(l2T)j) reduces to 

\I q \ < wTTl FT77 — VT ^ m 277-t 



91 " l-B (fei) S (k 2 ) ■ ■ ■ B (k m ) \ rr->o+ ' \B {k q )\ + & + - - • + £ m ) 2 + 1 

1 1 

x = 0, (21) 

\/(6 + 6 + --- + U 2 + i V^ + 1 

where £j £ [0, 1], for any j — 1, 2, . . . , m. Eq. ()21|) shows that 7 ? = for any g = 1, 2, . . . , m, 
and consequently ^j-y = for any E (lj). 

The fact that (3\ is maximal /minimal has played a key role in the arguments above. 
These arguments fail to hold if fii is not an extremum, since then the argument of 
5 (B (ki) x 2 + B (k 2 ) X4 + ■ ■ ■ + B (k m ) x 2m ) would have other zeros, besides the obvious 
x 2 = x^ = ■ ■ • = x 2m = 0. 

In the last step, we reveal the connection between the BEF for a general iV-level LZ- 



system and the exactly solvable bow-tie model J], in which all levels interact with only one 
special level (SL). Indeed, since S\\ does not depend on E if (3\ is maximal/minimal, 
we can safely set E (lj) = for any lj = 2, 3, ... , N. The form for Sn (Eq.flBJ) with all 
E (lj) = is exactly what one obtains for a bow-tie model if the SL has slope /3i, and a 
perturbation expansion is being developed for it. 

Since the Brundobler-Elser conjecture for the bow-tie model was proven analytically by 
Ostrovsky and Nakamura in 0, we only cite their statement from page 6947 of j^J: "This 
hypothesis is confirmed within the present model." , and refer to their work for further details. 
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This ends the proof of the Brundobler-Elser formula for a general iV-level LZ-system. 
Before drawing the conclusions, it is worthwhile to mention that, preceding [ij, an alternative 
)roof of BEF, based on analytic continuation in imaginary time, was proposed by Shytov in 



.1 



III. CONCLUSIONS 

In summary, we revealed the deficiencies of the BEF proof proposed by Volkov and 
Ostrovsky in yj and Appendix of Q. We emphasized that using the important property of 
the slope (3\ being an extremum only to arrive at the sub-series (5.1) in their proof is not 
enough, regardless of the changes of variables they attempted, and this property still has to 
play a crucial role in evaluating (5.1). 

We have constructed a proof of BEF that corrects for this shortcoming. Indeed, in our 
proof, the fact that the slope (3\ is maximal/minimal had to be used for three different 
purposes: 1) to arrive at Eq.(5.1) as was shown in jlj ; 2) to prove that Sn does not depend 
on E 3) to employ the fact that the Brundobler-Elser conjecture for maximum/minimum 
slope was proved for the particular case of the bow-tie model 

Note added. -After this work was submitted to publication we received a Reply by Volkov 



and Ostrovsky |8j in which these aut 



lors now verify, by direct calculation, the important 



point on which they missed before in l| and Appendix of j^J , and which we have advocated 
throughout our Comment, namely that in order for the BEF to hold the important property 
of Pi being an extreme slope {ES) must be used manifestly in the evaluation of series (5.1). 

In their Reply V-0 abandon the old changes of variables and use this time the same well- 
known change of variables (fTT]) (see, e.g., textbook |7[, page 114) as we have employed. We 
note that this change of variables, while previously used by V-0 in a different context (e.g., 
in the proof of the no-go theorem jjH, was never employed in their attempts of calculating 
series (5.1) in [l] and Appendix of|2|. The point at issue in our Comment is, however, not 
about what changes of variables should be or not be used in the proof of BEF. Mathematical 
inconsistencies aside, it was definitely not the V-O's choice of integration variables (in jlj and 
Appendix of 2]) that prompted our Comment, but their total neglect of the very important 
ES'-property of f3\ in evaluating series (5.1). 
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Though only the evaluation of 4-th order terms is presented in [8j] , our claim that one can 
set £j = only if (3\ is an ES is reconfirmed, together with the fact that this E'S'-property 
must still be used even after setting Sj = in order to prove the Brundobler-Elser conjecture. 
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